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In this article we compute the Casimir force between two finite- width mirrors at finite temperature, 
working in a simplified model in 1 + I dimensions. The mirrors, considered as dissipative media, 
are modeled by a continuous set of harmonic oscillators which in turn are coupled to an external 
environment at thermal equilibrium. The calculation of the Casimir force is performed in the 
framework of the theory of quantum open systems. It is shown that the Casimir interaction has 

£Nj , two different contributions: the usual radiation pressure from vacuum, which is obtained for ideal 

mirrors without dissipation or losses, and a Langevin force associated with the noise induced by the 
interaction between dielectric atoms in the slabs and the thermal bath. Both contributions to the 
Casimir force are needed in order to reproduce the analogous of Lifshitz formula in 1 + 1 dimensions. 

^H . We also discuss the relation between the electromagnetic properties of the mirrors and the spectral 

density of the environment. 
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I. INTRODUCTION 

Given the precision that has been recently achieved in the measurement of the Casimir forces [l| , the use of realistic 
models for the description of the media that constitute the mirrors is an unavoidable step for the improvement of 
Casimir energy calculations, which is needed for comparison with the experimental data. Moreover, from a conceptual 
point of view, the theoretical calculations for mirrors with general electromagnetic properties, including absorption, 
is not a completely settled issue @-|4| ■ Since dissipative effects imply the possibility of energy interchanges between 
different parts of the full system (mirrors, vacuum field and environment), the theory of quantum open systems [5| is 
the natural approach to clarify the role of dissipation in Casimir physics. Indeed, in this framework, dissipation and 
noise appears in the effective theory of the relevant degrees of freedom (the electromagnetic field) after integration of 
the matter and other environmental degrees of freedom. 

Dielectric slabs arc in general nonlinear, inhomogeneous, dispersive and also dissipative media. These aspects 
turn difficult the quantization of a field when all of them have to be taken into account simultaneously. There are 
different approaches to address this problem. On the one hand, one can use a phcnomcnological description based 
on the macroscopic electromagnetic properties of the materials. The quantization can be performed starting from 
the macroscopic Maxwell equations, and including noise terms to account for absorption [6(. In this approach a 
canonical quantization scheme is not possible, unless one couples the electromagnetic field to a reservoir (see [3[), 
following the standard route to include dissipation in simple quantum mechanical systems. Another possibility is 
to establish a first-principles model in which the slabs are described through their microscopic degrees of freedom, 
which are coupled to the electromagnetic field. In this kind of models, losses are also incorporated by considering a 
thermal bath, to allow for the possibility of absorption of light. There is a large body of literature on the quantization 
of the electromagnetic field in dielectrics. Regarding microscopic models, the fully canonical quantization of the 
electromagnetic field in dispersive and lossy dielectrics has been performed by Huttner and Barnett (HB) [7|. In 
the HB model, the electromagnetic field is coupled to matter (the polarization field), and the matter is coupled to 
a reservoir that is included into the model to describe the losses. In the context of the theory of quantum open 
systems, one can think the HB model as a composite system in which the relevant degrees of freedom belong to 
two subsystems (the electromagnetic field and the matter), and the matter degrees of freedom are in turn coupled 
to an environment (the thermal reservoir). The indirect coupling between the electromagnetic field and the thermal 
reservoir is responsible for the losses. As we will comment below, this will be our starting point to compute the 
Casimir force between absorbing media. 

Regarding the Casimir force, the celebrated Lifshitz formula Q describes the forces between dielectrics in terms 
of their macroscopic electromagnetic properties. The original derivation of this very general formula is based on 
a macroscopic approach, starting from stochastic Maxwell equations and using thermodynamical properties for the 
stochastic fields. As pointed out in several papers, the connection between this approach and an approach based on 
a fully quantized model is not completely clear. Moreover, some doubts have been raised about the applicability of 
the Lifshitz formula to lossy dielectrics 043- 

The first calculation of the Casimir force between two absorbing slabs using a microscopic approach is, to our 
knowledge, due to Kupiszcwska [9J, who modeled dielectric atoms as a set of harmonic oscillators coupled to an 
environment with T = , in which the atoms can dissipate energy. In that work, a scalar field in I+I dimensions 
was considered, and all the environmental effects were described through a dissipative constant and a Langevin 
force. In the context of quantum open systems, this is tantamount to consider an Ohmic environment. The force 
between slabs was then obtained in terms of the reflection coefficients associated to the slabs, which are described by a 
complex dielectric function. This result was rederived using a Green-function method for quantizing the macroscopic 
field in absorbing systems in ID in conjuction with scattering matrix approach [l0|. This was also extended to two 
identical absorbing superlattices [11] . Esquivel- Sir vent et al demonstrated an alternative Green- function approach that 
makes the quantization of the field within the slabs unnecessary and calculate the Casimir force in an asymmetric 
configuration [12| which was earlier considered only in the lossless case [ijjj]. In Ref.fl4J. the Casimir force was 
calculated in a lossless dispersive layer of an otherwise absorbing multilayer by using the macroscopic field operators 
considered by [15|] ■ In the series of works 0, Ha, LLZ| i Rosa et al considered the evaluation of the electromagnetic 
energy density in the presence of an absorbing and dissipative dielectric, using the HB model for T = and constant 
dissipation. In particular, in the recent Ref . [17| . they obtained the force density associated to spatial variations of the 
permittivity from which, in principle, one could obtain the Casimir force between slabs. 

In this paper we will follow a program similar to that of Ref. [9( , generalizing it by considering a general and well 
defined quantum open system. We will work with a simplified model analogous to the one of HB, assuming that the 
dielectric atoms in the slabs are quantum Brownian particles, and that they are subjected to fluctuations (noise) and 
dissipation, due to the coupling to an external thermal environment. We will keep generality in the type of spectral 
density to specify the bath to which the atoms are coupled, generalizing the constant dissipation model used in Ref. [9[ . 
Indeed, after integration of the environmental degrees of freedom, it will be possible to obtain the dissipation and 



noise kernels that modify the unitary equation of motion of the dielectric atoms. As we will see, general non-Ohmic 
environments do not provide constant dissipation coefficients in the equation of motion of the Brownian particles, even 
at high temperatures. Moreover, the spectral density of the environment determines the electromagnetic properties 
of the mirrors and therefore have a direct influence on the Casimir force. 

In addition to the conceptual issues described above, there are additional motivations to consider detailed micro- 
scopic models of the Casimir force, in particular, the controversy about its temperature dependence. Assuming simple 
phenomcnological descriptions of the materials based on the plasma or Drudc models, the theoretical predictions for 
the Casimir force are different due to the contribution (or not) of the TE zero mode [18| . At small distances a such 
that aT < 1, the differences are not too large, and the experimental results by Decca et al [19| seem to be well 
described by the plasma model. However, at large distances aT ^> 1, the theoretical predictions differ by a factor 2. 
The Casimir force at such large distances have been recently measured [20|, and the results are compatible with the 
Drude model after taking into account the interactions due to patch potentials on the surfaces of the conductors (some 
authors disagree with the evaluation of the effects of the patch potentials, see |2l[). In any case, these controversies 
show that more detailed microscopic models are necessary to clarify the situation. For example, considering that the 
slabs contain classical or quantum non relativistic charges interacting via the static Coulomb potential, the result for 
the large distance limit agrees with that of the Drudc model |22| . 

Another motivation for considering the Casimir forces in the framework of open quantum systems is the possibility 
of analyzing non equilibrium effects, such as the Casimir force between objects at different temperatures |23( and the 



power of heat transfer between them 24| , including the time dependent evolution until reaching a stationary situation. 
The paper is organized as follow: in next Section we shall present the model, the Heisenberg equations of motion 
for the different operators, and the vacuum and Langevin contributions to the field operator. In Section III we study 
much in deep the relationship between the microscopic model and the macroscopic electromagnetic properties of the 
mirrors. Section IV is dedicated to the evaluation of the Casimir force. After adding the vacuum and Langevin 
contributions, we show that the total force is given by a Lifshitz-likc formula, where the reflection coefficients of the 
slabs depend on the properties of the atoms and the environment considered in the model. In Section V we comment 
on the relation of the quantum open systems approach developed in this paper and the Euclidean computation of the 
Casimir force. We summarize our findings in Section VI. The Appendices contain some details of the calculations. 

II. THE MODEL 
A. The Lagrangian Density 

With the aim of including effects of dissipation and noise in the calculation of Casimir force, we will use the theory 
of open quantum systems, having in mind the paradigmatic example of the quantum Brownian motion (QBM) [5[. 

The model consists of a system composed of two parts: a masslcss scalar field and dielectric slabs which, in turn, are 
described by their internal degrees of freedom (a set of harmonic oscillators). Both sub-systems conform a composite 
system which is coupled to a second set of harmonic oscillators, that plays the role of an external environment or 
thermal bath. For simplicity we will work in I + 1 dimensions. In our toy model the massless field represents the 
electromagnetic field, and the first set of harmonic oscillators directly coupled to the scalar field represents the atoms 
in the slab. 

Considering the usual interaction term between the electromagnetic field and the ordinary matter, the coupling 
between the field and the atoms in the slab will be taken as a current-type one, where the field couples to the velocity 
of the atoms. The coupling constant for this interaction is the electric charge e. We will also assume that there is no 
direct coupling between the field and the thermal bath. The Lagrangian density is therefore given by: 

£ = £</, + Cs + C^-s + Cb + £-s-B 

= -d^d^fi + Anr) ( -mr 2 (x;t) mujlr 2 (x]t) ] + Airrie(j}(x;t)f(x; t) 

+ 4tt?7^ f -m n ql(x;t) - -m n ulql(x;t) J - 47T7?^ X n q n (x; t)r{x\ t), (1) 

n ^ ^ n 

where we have stressed the fact that r and q n are also functions of the position, i.e. each atom interacts with a thermal 
bath placed at the same position. We have denoted by r\ the density of atoms in both slabs. The constants A„ are 
the coupling constants between the atoms and the bath oscillators. It is implicitly understood that Eq.([T]) represents 
the Lagrangian density inside the plates, while outside the plates the Lagrangian is given by the free field one. The 
configuration of the slabs, of thickness d and separated by a distance a is shown in Fig.l. 
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FIG. 1. (Color online). Plates configuration denning the five regions in space. The plates are formed by oscillators which are 
in contact with thermal environments. All the system is in thermal equilibrium at a temperature T. In this work we consider 
a 1+1 version of this configuration. 

The quantization of the theory is straightforward. It should be noted that the full Hilbert space of the model H, 
where the quantization is performed, is not only the field Hilbert space H^ (as is considered in others works where 
the field is the only relevant degree of freedom) , but also includes the Hilbert spaces of the atoms H 5 and the bath 
oscillators Hb, in such a way that H = H<p <8> Hs <8> Hb- We will assume, as frequently done in the context of QBM, 
that for t < the three parts of the systems are uncorrelated and not interacting. Interactions arc turned on at t = 0. 
Therefore, the initial conditions for the operators <p 1 r must be given in terms of operators acting in each part of the 
Hilbert space. The interactions will make that initial operators to become operators over the whole space H . The 
initial density matrix of the total system is of the form: 

ftO) = &(0) ® js s (0) <8> jfe. (2) 

Since we are interested in the steady state (t — > +oo), we will assume thermal equilibrium at temperature T = 1/(3 
between the three parts. Each density matrix in Eq.® will be taken as thermal type (we have set h = ks = c = 1). 

B. Heisenberg equations of motion 

Starting from the Lagrangian (TT]), it is easy to derive the Heisenberg equations of motion for the different operators. 
They are are given by: 



Pn = m n q n , 



(3) 



p n = -m n u„q„ + X„r, 



(4) 



p = mr + eq 



(5) 



p = -mu%r + ^ ^™*i 



(6) 



D(f> = Anr/er, 



(7) 



where the operators p and p n are the conjugate momentum operators associated to the operators f and q n respectively 
Substituting Eq.© into Eq.®, we get: 



m n q n + m n uj n q n - X n r = 0. (8) 

As usual in the context of QBM, we solve the equations for the operators q n taking r as a source, and replace the 
solutions into Eq.([6]). In this way, the microscopic degrees of freedom in the slabs satisfy a Langcvin-likc equation of 
the form 

p= — mu^r — m— I dvy{t — t)t(x;t) + F(x;t), (9) 

dt Jo 

where the damping kernel 7 and the stochastic force operator F are the same as the ones of QBM (see [5| for a general 
and complete view). They are given by 

9 r+00 j( u ,\ 
1 (t) = — duj-^cos(ujt), (10) 



F(t) = £ -t^= L-^K + e^bl) . (11) 

Here b n and b^ are the annihilation and creation operators associated to q n , and J(u>) is the spectral density that 
characterizes the environment. This function gives the number of oscillators in each frecuency for given values of the 
coupling constants A„: 

J ^) = J2^r 6 ^-^- ( 12 ) 

ZTli n UJ r i 

In order to obtain a true irreversible dynamics, we introduce a continuous distribution of bath modes, that replaces 
the spectral density by a smooth function of the frequency u of the bath modes. Different functions will describe 
different types of environments. Physically, the thermal bath has a finite number of oscillators in a given range of 
frequencies. Then, a cutoff function must be introduced, containing some characteristic frecuency scale A. Then, the 
spectral density takes the following form 

In the equation above, 70 is the relaxation constant of the environment, while / is the frequency cutoff function. The 
values of a classify the different types of environments: a = 1 corresponds to an ohmic environment (in which there 
is a dissipativc term proportional to the velocity in the Langevin equation for the Brownian particle), while a < 1 
and a > 1 describe subohmics and supraohmics environments respectively [5[. 

At equilibrium, the stochastic force operator in Eq. ([TTj) and the damping kernel 7 in Eq. ([TU|) are not independent. 
The statistical properties of the stochastic force operator are given by the dissipation and noise kernels 

/r- - i\ r- - i r + °° 

D(t-t') = i( F{t);F(t') j =i F(t);F(t') =-2/ dwJ(w)sin[w(t-t / )], (14) 



D 1 (t-t') = {{F(t);F(t')})=2j 3 °d W JHcoth(^)cos[ w (t-t')], (15) 

which arc the formal quantum open systems generalization of the relations employed in Ref. [9j for general environments 
and arbitrary temperature. Note that only the noise kernel D\ involves the environmental temperature T as a 
parameter. Considering Eq. (fTU)) and Eq.fTTJJ, it is easy to show that: 



^-Mt-s) = --D{t- S ), (16) 

at m 

which relates the damping kernel 7 to the statistical properties of the stochastic force operator F. 

All in all, the set of equations to solve now are Eqs.©, ((TJ), and ((9J. 

It is possible to obtain a formal solution for the operators r{x; t) by considering the field as a source for the 
equation. This solution generalizes the crude approximation made in Rcf. |9j for the evolution of the microscopic 
degrees of freedom in the mirrors. It is given by: 

r(x;t) = G 1 (t)r(x;0)+G 2 (t)r(x;0) + — f dsG 2 {t- s) (f{x;s) - e$(x;a)) , (17) 

where G\$ are the Green functions associated to the QBM equation that satisfy: 

Gi(0) = l , Gi(0)=0, (18) 

G 2 (0)=0 , G 2 (0) = 1, (19) 

for which, the Laplace transforms are given by 

G x {z) = - =? — — -, (20) 

G 2 (z) = - i — — -, (21) 

y ' z 2 + oj 2 + zj(z) 

where 7 is the Laplace transform of the damping kernel. Note that, given these conditions, one can prove that 
G x {t) = G 2 (t). 

Inserting this solution into Eq. ((7J), we obtain the following equation for the field operator 

□0 + _JI!Z£_ / Gift - tU(x; r)dr = 4wr]e \G 1 (t)r(x; 0) + G x (t)?(x; 0) + — / G x (t- t)F(x; tUt) , (22) 

m J Q V m Jo J 

subjected to the free field initial conditions: 

1 

> ' -~~ ikx I <>t —ike 



(x;0)= J dk[ — ) [a k e lkx +ale-^ x ), (23) 



i{x; 0) = fdk ( j- J 2 (-iu k a k e ikx + iu k a\e- ikx } , (24) 

where a k and a k are the annihilation and creation operators for the free field, and uj k = \k\. The boundary conditions 
are the continuity of the field and its spatial derivative at the interface points. 

We will compute the Casimir force from the xx-componcnt of the energy-momentum tensor 

f xx (x;t) = ±((d Q 0) 2 + (dJ) 2 y (25) 

The force is explicitly given by 

Fc = (T-*) - (f^>, (26) 



where the expectation values are taken on the regions outside the planes (Regions I or V) and between them (Region 
III), respectively, in a thermal equilibrium situation. 

For this calculation, we need the explicit solution for the field equation (|22"|) . Considering the properties of the 
fundamental solutions Gi,2, the first step is to take the Laplace transform of the equation in order to obtain 

2 (x; s) - s 2 ( 1 + ^^G 2 (s)) 4>(x- s) + s0(x; 0) + $(x; 0) + ^Gi(s)%; s) + ^Gi(s)p(x; 0) 



9a; 2 

+47r?7ef (ar; 0) (sG 1 (s) - l) = 0. (27) 

Since we are interested in the long-time behaviour (t — > +oo), we can omit terms containing the positions r and 
momenta p of the oscillators at t = 0. This assumption is well justified for t ^> I/70, the scale associated with the 
damping or relaxation time of the environment. Therefore, the equation for the field can be approximated by: 

^U(x; s) - a 2 (l + ^-G 2 {s)\ ${x; s) = - 8 $(x; 0) - $(x; 0) - ^Eq x (s)F(x; a). (28) 

This is the equation for the Laplace transform of the field, with the initial conditions and the Laplace transform of 
the stochastic Langevin force as sources. For simplicity, the spatial dependence of the matter terms was omitted, but 
it is important to remember that this expression is valid for the points inside the plates. Note also that, as there is 
no contribution from the atoms to the sources, the field operator <f> acts on Hs as an identity. 

We propose a solution of the form 4>{x;t) = <pv(x;t) + 4>l{x;£), where two contributions are distinguished: the 
vacuum contribution <J)v{x;t), which results from the modified field modes; and the Langevin contribution <^,(x;£), 
which depends linearly on the Langevin forces. Each contribution satisfies 

<f> v (x; s) - s 2 ( 1 + —G 2 (s) 4>v{x] s) = -s0(.t; 0) - <j>(x\ 0), (29) 



dx 2 



m 



^H{x- s) - s 2 (l + ^-G 2 {s)) M^ s) = -^Gi(s)%; a). (30) 



dx 2 ' \ m 

It is worth to note that the first equation presents only operators acting on H^,, thus the associated field contribution 
is an operator on that space. In the same way, the second equation depends only on operators acting on Hb, and 
therefore the Langevin contribution acts nontrivially only there. Taking all this into account the field operator reads 

4>(x;t) = 4> v (x;t)®I s ®I B +I < f,®l s ®'$L{x;t). (31) 

In summary, the atoms act like a bridge between the field and the thermal bath, making no contributions to the 
total field. The problem, then, has been completely separated into two parts, that will be computed in the next 
subsection. 

C. Vacuum and Langevin contributions 



We will now solve the two independent Eqs. (|29| and (|30|). For the vacuum contribution </>y, the solution for 
t — > +00 is assumed to have the form: 



4>v(x;t) 



+ OO 







d i& +3 » f " 



— — «->'«.) + «[<«• A*W , (32) 



where the first integral comprises the waves going from left to right, and the second one the waves going from right 
to left. The mode functions fk(x) satisfy the equation: 

^fk(x) + u 2 n 2 {u k )f k (x) - 0. (33) 



The refractive index n(uj k ) given by 

n 2 (u Jk ) = l + ^^G 2 (-iu Jk ) = l + - , 4 , (34) 

m luq - w£ - iu k "f(-iu k ) 

where Wp = nr,e is the plasma frecuency. 

It is worthy to note that Eq. (f3"3")l is of the same form that the equation for the modes in a nonabsorbing dielectric 
medium [26|, except that in this case the refractive index is frecuency-dependent. On the other hand, Eq. (|34[) depends 
on the Fourier transform of the damping kernel, which is associated to a general environment. There is no assumption 
neither about the spectral density of the environment nor on the value of the equilibrium temperature. In the next 
Section we will analyze in more detail the relation between the electromagnetic response of the medium and the 
spectral density. 

The Eq. (|33|) for the modes can be solved by imposing continuity of the mode functions f k at all interfaces. The 
calculation is long but straightforward and is described in Appendix A. 

In order to obtain the Langevin contribution, we must solve Eq. (|30p . Since the interaction begins at t = 0, the 
equations for the Laplace and Fourier transforms for this contribution are identical, and related by s = —ik. Thus, 
the long-time solution can be written as <pi,(x;i) — J_^ t^4>l{x; k)e~ lkt , where 4>l(x; k) satisfies 



dx 2 
in regions I, III and V; and, 



d 2 ~ 
—^(f>L(x;k)+k 2 ^L(x]k)^0, (35) 



9 2 -? / ,n ,o 2/n? / ,\ 47T7ye ik 

— M x;k) + k n W^(x;fc) = — ^ 2 _. fc ; , ( _. fc)) . 



5-20i(s; fc) + k 2 n 2 (k)cj, L (x; k) = '- ———F(x; k), (36) 



in regions II and IV. Here F(x; k) is the Fourier transform of the stochastic force operator. The explicit solution is 
presented in Appendix B. 

In Section IV we will use the vacuum and Langevin contributions to the field operator in order to obtain the 
Casimir force between slabs. Before doing that, we will describe in more detail the relation between the macroscopic 
electromagnetic properties of the slabs and the microscopic model. 

III. GENERALIZED PERMITTIVITY FROM THE QUANTUM OPEN SYSTEM 

With the aim of checking if our model is physically consistent, we will analyze the properties of the refraction index 
given in Eq. (|34[) . Considering that e(w) = n 2 (w) is the permittivity of the material plates, we have: 

e(u)-l = - 2 2 Ul . =p . . . (37) 

We can define the susceptibility kernel x( T ) f° r the model as [25( 

i r+oo 2 r+oo „-iut 

X(r) - 5- / (c(«) - lje-^dw = / / 2 . -, . . cL>. (38) 

In principle, this integral can be evaluated by contour integration. Inversely, the permittivity can be expressed in 
terms of x( T ) as: 

/+oo 
X (T)e^dr, (39) 

-oo 

which can be viewed as a representation of e(oS) in the complex w-plane. The permittivity is well defined when x( T ) 
is finite for all r and xi T ) -> as t -> ±oo, and its analytical properties can be studied directly from this expression. 
All properties of the permittivity and susceptibility functions are strongly dependent on the Laplace transform of 
the damping kernel ~f(ik), which in turn depends on the spectral density of the environment. After taking the Laplace 
tranform of Eq. (TITj)) , we obtain 



7(s) = ~f "d^TT^ (40) 

w m J u (s 2 +oj 2 ) v ; 

As already described, a physical spectral density must incorporate a cutoff function. One could use a sharp cutoff or, 
alternatively, choose a continuous cutoff function that approaches zero rapidly for frequencies greater than the cutoff 
frequency A, ensuring the convergence of the integral. 

The first alternative, although simpler, makes the function j(ik) not well-defined in the complex plane. The second 
alternative solves this problem and allows the use of the residue theorem to evaluate the integral. Inserting Eq. (TT3]) 
into Eg. (J4T)|). we obtain 

700 = 4^1 P^T-^Wfx)- (41) 

w ttA"" 1 J (s 2 +uj 2 )\AJ y ' 

In order to apply the residue theorem the integrand must be holomorphic on the superior complex half-plane, 
except in a finite number of points which are not on the real axis. Thus, different results are obtained considering 
distributions with or without poles on the superior half-plane. For an ohmic environment (a = 1) and no cutoff 
function it is easy to see that j(s) = 2j - 

In the case of distributions without poles (for example, a gaussian cutoff function), we have: 

j N p(-ik) = ^J(k)=j 1 (k), (42) 

where the subscript NP denotes the fact that the distribution has no poles. The resulting function is real and even in 
the variable k. 

On the other hand, the distributions usually considered in the literature have poles on ±iA (for example, a Lorentzian 
distribution). In these cases, for odd uj a (with a < 4 to mantain the convergence in Eq . (|4"Tj) ) . we get 

7p(-*) = ^AW+«(-l)^^(4) a J-k(A), (43) 

where the subscripts on J denote the location of the pole. Although the resulting function is complex, the second 
equality in Eq. (|42[) remains valid. 

Taking into account the above properties of the damping kernel, we now continue analyzing the properties of 
the permittivity and susceptiblility functions. As a particular example, in the Drudc model one has y{— iui) = 70- 
Therefore the denominator in Eq. (|38[) has two poles, both in the lower half w-planc. Thus, as expected from a physical 
point of view, the susceptibility kernel shows a causal behavior, since it vanishes for r < 0. The analyticity of e(w) in 
the upper half w-plane allows the use of Cauchy's theorem, resulting in the well-known Kramers-Kronig relations for 
the real and imaginary part of the permittivity function e(w). 

In our more general case, the physical properties of e(w) are determined by the function j(— iu>). This dissipation 
function is given by the theory of quantum open systems through Eqs. (T4"2"|) and (|4"3"|) . and depends on the chosen cutoff 
function. 

Let us first consider, for simplicity, the case in which the cutoff function has no poles, represented by the Eq.([4"2"j). 
For a given spectral density, the denominator in Eq. (|38p reads 



H=^-^ 2 -*27o^(^f '/QJ). (44) 






If we choose an ohmic environment (a = 1) and no cutoff function (which is equivalent to put /= 1), we reobtain 
the Drude model (if wo = 0) or the one-resonance model (when wo 7^ 0). In principle, we could consider other values 
of a while keeping / = 1. In this case, u° should be an odd function. For example, a = 3 gives an ill-defined 
pole configuration, since one of the poles lies on the upper half w-planc, breaking the analyticity of the integrand in 
Eq. (|38l) . and resulting in a non-causal susceptibility, which turns out to be unphysical. 

Therefore we see that for this supraohmic environment the use of a cutoff is unavoidable. We may use an analytical 
cutoff (like a gaussian function), or a lorentzian cutoff function. The first alternative leads to a denominator D^ p 



10 

whose zeroes cannot be obtained analytically. The second alternative, valid for a < 4 such that u> a is an odd function, 
leads to a denominator 



4 al M 



which for a — 3 gives 



(A 2 + uj 2 ){uj 2 - uj 2 ) + 2 70 A 3 - Q W 2 ((-l) i 
(A 2 + uj 2 ) 

Aw 2 — iu>lu - (270 + A)cj 2 + iuj 3 



DP(u) 



A — iui 



(45) 



(46) 



We denote the zeroes of D p ' as uii = u>oXi (with i = 1, 2, 3). The three roots turn out to be located in the lower half 
w-planc, which ensures the causality property. Also, one of the roots is purely imaginary (xi = —x* — —i\x\\) and 
the others two have the same (negative) imaginary part but opposite real parts (#3 = — x%). Thus, the susceptibility 
kernel reads 



,( 3 ) 



x£V) 



Up 
Ld 



(A-ajo|gi|) - Ua \xi\r 



(xi - x 2 )(xi +x%) 



2Re 



(A — 1UJ0X2) 

(X 2 ~ Xi){x 2 +£2) 



(r), 



(47) 



where it is clear that it is a causal real function and, due to the negativity of the imaginary part of the roots Xi , we 

have Xp ( T ) "^ f° r T ~ * +°° as is expected. We also have that Xp (0) = but %p (0) ^ 0, and therefore the 
asymptotic expression found in [2^] still remains valid, as well as the Kramers-Kronig relations. 

It is worth noting that an ohmic environment (a = 1) can be also studied with a cutoff function, obtaining similar 
results. 

All in all, we have shown that our model is physically consistent and generalizes previous results for the permittivity 
of absorbing media, including as particular case the Drude model. Plasma-like models do not contain dissipation and 
can be obtained by taking 70 = 0, that corresponds to no coupling between the system and the bath. 



IV. CASIMIR FORCE 

A. The Energy-Momentum Tensor and the Different Contributions to the Casimir Force 

Once determined the two contributions to the field, we proceed to compute the Casimir force between the plates 
as given by Eq. (|26[) . For this purpose, it is necessary to compute the expectation value of T xx which is given by Eq. 

Considering the vacuum and Langevin contributions according to Eq. (l3ip . we have 



T, 






(48) 



It is worth to remark that there are cross terms which act over two parts of the total Hilbert space. 

As we are interested in the steady state of the system, which is assumed at thermal equilibrium, each part of the 
total density matrix is represented by a thermal-type density matrix. On the other hand, both field contributions 4>v 
and 4>l are linear on the annihilation and creation operators of their respective parts of the total Hilbert space. Thus, 
the cross terms do not contribute to the force in the case of thermal equilibrium. Then, the problem is reduced to 
compute the expectation values over thermal states of the operators T^ x and T£ x , i.e., 



j=(T xx )= Tr^ [p+fn + Tr B (p B T^ I = f v + f L . 



(49) 



Thus, the Casimir force also have two contributions: 



Fc=h- fin = (fX + tf) - (fYn + ff 11) =F% + F£. 



(50) 
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B. Vacuum Casimir Force 

For the vacuum contribution, T xx is quadratic in the annihilation and creation operators. Thus, in order compute 
the expectation value over the thermal state, we need to evaluate the expectation values of the products of the 
annihilation and creation operators. These arc given by the known expressions: 



{a k a k >) = (a k al,) = 0, 



(51) 



a k a{,} = 6{k-k')(l + N{cj k )), 



(52) 



(a{a k ^ = S{k-k')N(u} k ), 



(53) 



where N(u) k ) = efil } k _ 1 . 

Taking into account Eq.(|32j). we have, in region I 



fr(x)=Tr^[p 4> T^) = - 



1 



o 1 



d ±r^(^\ 



^coth^n. fc |/<(,)P + -l 



df, 



d.v 



(54) 



which is identical to the expression for a non-absorbing medium except that, in this case, there is a thermal factor 
coth(/3a>fc/2) related to the temperature of the field. 

Using the solutions for the modes functions f k in regions I and III (see Appendix A), the vacuum Casimir force is 
given by 



*c —h - tin - 



\ \T £ fccoth (ir) t 1 + |i?fc|2 + lTk]2 2 ( |Cfc|2 + lDk ^ 



(55) 



The coefficients R k , T k , C' k , and D k are explicitly given in Appendix A. It is worth noting the appearance of a 
thermal global factor in the last expression, which comes from the field's thermal state at temperature T, based on 
the equilibrium assumption. 



C. Langevin Casimir Force 



For computing the Langevin contribution to the force, it is necessary to know the expectation value, over the bath's 
thermal state, of the force operator, where binary products are evaluated at different frequencies. 

For any time-dependent hcrmitian operator, the expectation value evaluated at different times, corresponds to 
the correlation function of the operator. This matches with one-half of the anti-commutator expectation value at 
different times. Thus, making Fourier transforms over both times, we can compute the desired products of the Fourier 
transform of the force operator, at different frequencies. 

For the case of thermal equilibrium, the anti-commutator expectation value at different times of the force operator 
is provided by the QBM theory. One can show that it matches the noise kernel D\{t — t') of Eq. ([T5)) . Thus, we obtain 



{%);%')}) = JK)coth (^f) 



V 2 



(56) 



Considering that in our case the stochastic force operator depends on the position where the atom is located, we 
finally have 



,=S(x- x')^^- coth (^1) S(k - 



(57) 



where we have included the atom density 77 due to dimensional issues. As it can be seen, the frequency spectrum is 
not flat. 

Taking all this into account, T xx in regions I and III are given by: 
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■*■ xx \ x it) 



Bl / ^-kk')W 1 (k)W 1 (k')e-^ k+k ^ +t \ 



2n 



2tt 



(58) 



hLJII 



(x;t) 



+°° dk 

.r>o 27T 



dk' /^~~ ^~~ ^~~ -~— 

{~kk')e~ l{k+k )l [W 2 {k)W 2 {k')e l{k+k > + W 3 (k)W 3 (k')e- l{k+k )x 

2tt V 



(59) 



The coefficients Wi(k) in these equations arc given in Appendix B, and arc linear functions of the Fourier transform 
of the stochastic force operator. Taking into account the explicit expressions in Appendix B and Eq. (|57[) . the desired 
expectation values are: 



(Wi(fc)Wi(fc')) = I2»( fc )l 



,27T J(Wfc) 



coth 



f /3wk 



5{k + k') 



m /c 2 7i(fc) V, 2 

m (1 - e- 22 ) (|£| 2 e 22 + |r„ + re l2ka \ 2 + e 22 |l + rr n e l2ka \ 2 + \t\ 2 \r n \ 2 ) 

+ 2n 2 Re (\t\ 2 i (e~ lZl - l) r* + * (l - e 121 ) (l + r*r* n e~ l2ka ) (r„ + re l2fea )) 



(60) 



coth 



ni (1 - e~ 22 ) (1 + |r„| 2 e- 22 ) + 2n 2 e~ Z2 Re (i (l - e iZl ) r„) 



• = (t? 3 (fc)^ 3 (fc')) = l^^l 2 ^^!^ (! 



5(jfc + A-') 



(61) 



where n = rii + in 2l i.e., m = Re(n) and ri2 = Im(n), ^/(ik) = 7i(fc) + i^ 2 (k), z\ = 2kri\d and 2:2 = 2kn 2 d. The 
explicit expressions for the coefficients r, r„ and t can be found in Appendix A. 
Therefore, the Langevin contribution to the force in regions I and III are given by 



jf=Tr B [p B T 



L XX 



? »)l 2 -*T 

Air m 7i(fc) 



coth 



(¥) 



x ni (1 - e- 22 ) (|t|V 2 + |r„ + re l2fcQ | 2 + e 22 |l + rr n e aka \ 2 + \t\ 2 \r n \ 2 ) 
+ 2n 2 Re (\t\ 2 i (e^ - l) r* + * (l - e< 21 ) (l + r*<e" j2fca ) (r„ + re l2fca )) 



(62) 



fin - 2> B p B 3^" 



" dk 

2^ 



frer 4lJK) 



coth ( — — 
m 71(A) \ 2 



H (1 - e- 22 ) (1 + |r„| 2 e" 22 ) + 2n 2 e" 22 .Re (* (l - e 121 ) r„) 



(63) 



Taking advantage that the integration is over all the values of A;, the fact that the change k <R- — k is equivalent to 
complex conjugation, and the second equality of Eq. (|42|). we obtain 



T/-y o 27T2 |n|2 |l_ r 2 e i2fca|2 COtJl ^ ? jl 1 



D i2fca|2 , z 2 



e 2 1 + rr„ e 



i2fca|2 



+ l*r>»f 



(64) 



/• +00 dfcfc|n + l| 2 16|t| 2 m 



II 



2tt2 \n\ 2 \l-r 2 e 



1„i1ka\1 



Bk 

2 



i + mi + kre 



2 c -^2^ 



(65) 



Note that the presence of the thermal factor coth ( 4p ) is in agreement with the null temperature result obtained 



in other works for an ohmic environment |9| since when T — > 0, coth 
Finally, the Langevin contribution for the Casimir force is: 



^C — 7/ - )iii — 



+ °°dfcfc |n+l| 2 S\t\ 2 m (f3k 

2tt2 \n\ 2 \\- r 2 e i2kaV2 COi I 2 



x (1 - e- 22 ) (|i| 2 + \r n + re l2ka \ 2 e^ + |1 + rr n e l2ka \ 2 + \t\ 2 \r n \ 2 e - z * - 2 (l + \r\ 2 ) (l + |r„| 2 e - 22 )) . (66) 
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It should be noted that here also appears a global thermal factor, as in the vacuum case, but this comes from 
the bath's temperature while in the vacuum case comes from the field's equilibrium temperature. Note also that Fq 
vanishes when there is no coupling to an environment, since in this case the refraction index is real and therefore 

22=0. 

D. Total Casimir Force 

The total Casimir force is determined from the expressions Eq. (|50p . ([55)1 and (p5|) . The resulting force can be 
written in a very compact form. On the one hand, it can be proved that the total free energy in region I (outside the 
plates) coincides with that for free field at temperature T i.e., 

''=«'+'' =r°^» th (??)> < 67 > 

which is expectable due to translational invariance outside the plates and our assumption of thermal equilibrium. On 
the other hand, in region III we have 



X dk (Pk\ (l-|r| 4 ) 

2tt 



hn = fin + tf h =J q 



Therefore, the total Casimir force is finally written as: 

dk 

— .„„ -_^ 2 J \" |l_ r 2 el 2fca|2 



+ < "g fcco th(^) 



or equivalcntly 



1 



TV 



/3k\ r 2 e 



■1 



2A2ka 



F c [a] = --Re rffcfccoth £- - , ... . (70) 



This expression is formally identical to the case of dissipative material (ohmic environment) at zero temperature 
found in previous works [9|, l26|. However, it contains two generalizations: on the one hand, temperature has been 
included in the formalism in a natural way by means of the open quantum system theory. On the other hand, the 
equation is valid for a general environment: the refraction index is in general complex, and dependent on the function 
j(ik), which comes from the interactions at the microscopic level. The associated permittivity depends on the type 
of environment, and reproduces known results (as the Drude model) as particular examples. 

E. Convergence and Limits 

In this Section we will study some properties and limits of our final result Eq.(|70p. Let us first study the convergence 
of this expression. In general, the Casimir force calculations involve several rcgularization methods to achieve a finite 
result. A usual approach is to introduce a high frequency cutoff, in order to take into account that real materials 
are transparent at high frequencies. This characteristic is already incorporated in our model. Indeed, the complex 
refraction index n includes all the environment properties which produce dissipation and noise. For large values of k, 
taking into account Eq.([34|) one can check that 

"i^ 1 -— —73. ( 71 ) 



27T7 ? e 2 7i(fc) 

n 2 -> ro-. (72) 

m k^ 

Then, in the same limit the reflection coefficient r behaves as 

^—4(1-^), (73) 
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and in consequence the integrand of Eq. ([70|) is 0(k~ 3 ) for large values of k. Thus the convergence is ensured when 
k — > +00, regardless the temperature and the type of environment considered. 

On the other hand, Eq.([70|) contains as particular cases some known results. The non-absorbing medium case 
can be easily obtained by setting the relaxation constant 70 = in all the expressions. This makes the refraction 
index real, which cancel the Langevin contribution in Eq. ([66p since the factor 1 — e~ 22 — > 0. Thus, only the vacuum 
contribution in Eq. ()55[) survives but with a real refraction index [26J . 

Another important limit is the well-known Lifshitz formula. In the original work [8[ , Lifshitz considered semispaccs 
separated by a finite distance. Therefore one should take the limit of large thickness (d — > +00), in which r should 
be replaced by r„. After a rotation to the imaginary- frequency axis Eq. (|70[) becomes, in the T = case 



1 f + °° r 2 (is)e~ 2sa 

which is of the form of Lifshitz formula for a scalar field in 1 + 1 dimensions. In the case T ^ 0, wc must take into 

account that coth I ^- ) has poles on the imaginary axis at the Matsubara frequencies 2nif3j = i£j,j = 0,1,2,.... 

Therefore, the integration path can be rotated to the imaginary axis, but must be deformed to avoid the poles. This 
is a standard procedure, that converts the integral over frequencies into the Matsubara sum 

■^ r 2 (if-)e^ 2 ^ a 
Fc[a] = 2T^£ j r "f. )e (75) 

which is the standard expression for Lifshitz formula at T 7^ 0. 

It is interesting to remark that in our simplified 1 + 1 model there is no discontinuity in the transition between the 
Drudc and plasma models. The Drude model is recovered, in the ohmic environment case, by setting ujq = (i.e. free 
particles instead of harmonic oscillators for modelling the dielectric atoms) and, once this limit is taken, the plasma 
model corresponds to the particular case 70 = (no coupling to the environment). The reflection coefficient r n — > 1 
in the zero frequency limit, for any value of wp and 70, even setting 70 = from the beginning. In the case of absence 
of coupling to an environment, of course one must assume thermal equilibrium. Analogies between thermodynamics 
of a free Brownian particle and that of an electromagnetic field between two mirrors of finite conductivity have been 
studied in Rcf. [2Z|. 

V. CONNECTION WITH THE EUCLIDEAN FORMALISM 

Given that we are assuming thermal equilibrium between the different parts of the system, the results presented in 
this paper for the Casimir force could be derived following a functional approach in Euclidean space, as described for 
instance in Refs. |28|,l29J]. We mention briefly the relation between both approaches. 

As is well known, in 1 + 1 dimensions the free energy £ for a quantum system in thermal equilibrium at temperature 
T can be computed as 

f=-Tlog|£) (76) 

where Z(a) is the partition function when the plates are separated by a distance a. The partition function can be 
represented by the functional integral 

Z = fv(f>VrVq n e~ SE , (77) 

where Se is the Euclidean action for the full system. The integration is performed by imposing periodic boundary 
conditions on the temporal coordinate. 

After integrating the matter and bath degrees of freedom it is possible to find an effective action for the scalar field 
of the form 

S cS = J (Px-d^d^ + J d 2 x J d 2 x'V(x, x')<f){x)<f){x') , (78) 

and the vacuum persistence amplitude reads 

Z = \V<f>e- s * t{ . (79) 
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The effective action for the scalar field is quadratic because of the linear coupling we arc choosing for the interaction 
between the vacuum held and the matter degrees of freedom. The potential V(x,x') is different from zero only inside 
the plates, and as in our model the field 4>(x, t) interacts only with the atom at x, it can be shown that the potential 
is of the form V(x, t, x',t) = 8{x — x')X(t — £'). The function \(t — t') encodes the information about the interaction 
of the vacuum field with the matter degrees of freedom, and also of the influence of the thermal bath. Its Fourier 
transform is related to the reflection coefficient of the slab. 

Formally, the vacuum persistence amplitude is given by the functional determinant 

Z = (det[-n + V])~* ■ (80) 

An explicit evaluation of this determinant leads to Lifshitz formula [29( . So when considering thermal equilibrium, 
one has an alternative route to the evaluation of the Casimir force, even when the field is coupled to other degrees of 
freedom. However, this Euclidean functional approach would not be adequate to compute the force for more other 
initial states or, in general, in noncquilibrium situations. In that cases, the use of the theory of quantum open systems 
as described in this paper is unavoidable. 

VI. CONCLUSIONS 

In this paper we have presented a derivation of the Casimir force between two absorbing slabs in the framework of 
the theory of quantum open systems. We worked with a simplified model of a scalar vacuum field in 1 + 1 dimensions. 
In order to describe the interaction of the vacuum field with the mirrors, we considered a model analogous to the HB 
model for QED, where the matter degrees of freedom are described by a continuous set of harmonic oscillators, which 
are coupled not only to the vacuum field but also to a thermal bath that accounts for dissipative effects. 

Following a standard procedure in the theory of quantum open systems, we showed that the field operator satisfies 
the modified Klein Gordon equation (|22|) . This is a nonlocal Langevin equation, which describes the interaction of the 
vacuum field with the matter degrees of freedom and the effects of the thermal bath on its dynamics (the environment 
is indirectly coupled to the quantum field through the matter). This equation is similar to the one that describes a 
Brownian particle coupled to an environment. Both the noise and dissipation are determined by the properties of the 
environment. 

The field operator that solves this "Klein Gordon- Langevin" equation can be written as the sum of two terms, a 
vacuum contribution and a Langevin contribution. The same happens with the associated energy-momentum tensor, 
and therefore we have a similar decomposition for the Casimir force between slabs. The final result for the total 
Casimir force is equivalent to a 1 + 1-version of the Lifshitz formula, expressed in terms of the reflection coefficients 
associated to the slabs. Therefore, we have presented, in this simplified model, a first-principles derivation of Lifshitz 
formula in the framework of quantum open systems. 

The present work is closely related to Ref. [3], that has been improved and generalized in several directions. Indeed, 
that work assumes the simplest forms for noise and dissipation (ohmic environment and constant dissipation) , without 
specifying the properties of the environment. Moreover it is doubtful whether a general non-ohmic environment can 
produce such effects at T = 0. Here we worked at T ^ 0, and considered very general environments. We also linked 
the properties of the environment with the macroscopic electromagnetic properties of the mirrors. There is also a 
close relation with the recent work [17J , where the authors computed the force density associated to spatial variations 
of the permittivity. As compared with ours, in this reference the authors considered the more realistic case of 3 + 1 
electromagnetic field, but only in the particular case of T = and constant dissipation. Moreover, they did not 
consider the presence of boundaries as we did here, which allowed us to compute explicitly the Casimir force between 
slabs and to obtain Lifshitz formula. 

In order to apply the quantum open systems approach to a realistic calculation of the Casimir force, we should 
generalize our results to a 3 + 1 model with the electromagnetic field. Although technically more complex, we do 
not expect conceptual complications in doing so. Regarding the long standing controversy about the temperature 
corrections to the Casimir force, a crucial point is the behavior of the quantity 

limC 2 [c(»C)-l], 

which vanishes for the Drude model and is different from zero for the plasma model, producing in the last case an 
additional contribution to the force coming from the TE zero mode. In the kind of microscopic models considered 
here, the TE zero mode is suppressed as long as ujq -£ 0, as can be seen from Eq. (|3"7|) . 

Finally, in the 3 + 1 dimensional models one could consider more general initial states and/or non-equilibrium 
situations. We hope to address this issue in a forthcoming publication. 
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Appendix A: Boundary conditions and solutions for the vacuum contribution 

In this Appendix we present the explicit form of the solutions of Eq. p3[) . En each region we have: 

fi(x) = e lkx + R k e- ikx , (Al) 

fl*(x) = A k e lknx + B k e- lknx , (A2) 

f I k II {x) = C k e: kx +D k e- lkx , (A3) 

fV( x ) = E k e tknx + F k e~ lknx , (A4) 

f k v (x)=T k e lkx . (A5) 

The different coefficients can be obtained by imposing continuity of the mode functions and their derivatives at the 
interfaces. They are given by: 

/ 4.2 p i2ka \ 

i? fc = e-^(j-+ (1 _ rV2fca) ] ; (A6) 

t(\ -4- r rP i2ka \ 

A k = e- ife *e ifc "* U , rL , (A7) 

t n {l-r 2 e t2ka ) v ' 



ik%„ikn% K r n + r£ ) 

t n (l - r 2 e i2ka ) : 



B k = e-^^ :^- r ';> ( Ag ) 



t 

(\ — r 2 R i2ka\ ' 



C « = X ^i2ka, > ( A9 ) 



tre ika 

D * = {l _ r 2^2kay ( A1Q ) 



TT, — . ik§ ikd -ikn§ -iknd (Alii 

*"t n (l-r¥ 2 *») ' [ ' 

Fk = , l r " ■■». e ^ e lkd e tkn ^e lknd , (A12) 
t n (l — r z e lZka ) 



Tk ~ ft r 2 p i2kaY ( A13 ) 



(l __ j.2 c,i2ka\ '■ 

for fc > (while for k < the order of the solutions must be reversed and the refractive index and the coefficients should 

/ i2knd -i\ . iknd ikd 

be conjugated), where r = ,"_ 2 ;;,,„,, and £ = , Vlyi if 2 e i2knd\ are the reflection and transmission coefficients for 

one plate, while r n = ^^r an d t n = -^V are the ones for an interface. The coefficients R k and T k can be interpreted 
as the reflection and transmission coefficients of the two plates configuration. However, it should be noted that, due 
to the presence of absorption, \r\ 2 + \t\ 2 7^ 1 and |-Rfc| 2 + \T k \ 2 7^ 1 in this case. 
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Appendix B: Boundary conditions and solutions for the Langevin contribution 

In this Appendix we solve the equations (|35|) and (|36|) . Due to the presence of a source in regions II and IV (see 
Eq. (|55|) ). the solutions will have to parts: one associated to the homogeneous equation and other related directly to 
the source. Therefore, the solutions are: 

ft L (x;k) = W 1 (k)e- lkx , (Bl) 

^{x\ k) = t/i(fc)e 4fc " :r + U 2 (k)e~ lknx + A t {x\ k)e tknx + A 2 (x; k)e~ iknx , (B2) 

fi L "(x; k) = W 2 (k)e lkx + W 3 (k)e~ tkx , (B3) 

ft/ix; k) = V!(k)e lknx + V 2 {k)e- %knx + B x (x; k)e tknx + B 2 (x; k) e - tknx , (B4) 

^(x;k) = W i (k)e ikx , (B5) 

with, 

A 1 (x;k) = -!- I G(x'-k)e- lknx 'dx', (B6) 

Aa(x]k) = ~ I G(x';k)e iknx 'dx', (B7) 

2n J_ d _a 

B 1 (x]k) = ^- I G(x';k)e- tknx 'dx', (B8) 

B 2 (x;k) = - — I G{x';k)e- tknx 'dx', (B9) 

where, for simplicity, we write G(x; k) = 4^__gg__. 

The coefficients Wi(k), Ui(k) and Vi(k) are obtained by means of the appropriate boundary conditions. Thus, they 
are given by: 

Wi(k) = W{k)e tknd e-^ a+d ^ (k (1 + rr n e l2ka ) +L(r n + re l2ka ) + Mte lkd e tk ^ a ~ nd ^ + Ntr n e tkd e lk( - a+nd A , (BIO) 

Wa(fc) = OJ(fc) (Kr n e l2knd + L + Mre lka + Nrr n e lk( - a+2nd] ) , (Bll) 

W 3 (k) = W(k) (krr n e tk{a+2nd) + Lre tka + M + Nr n e l2knd ) , (B12) 

fj x = !Jl e iKn+m%+d)w u (B13) 

U 2 = —eW-^M+QWi, (B14) 
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where, 



Vi = — e tt(1_T,) *(W 2 +r n e- lka W 3 ), (B15) 

V 2 = ^-e 1 ^ 1 - 1 ^ {r n e lka W 2 + W 3 ), (B16) 



W 4 = e ife ("- 1 ) d (W 2 + r n e- %ka W 3 + - -N), (B17) 

(n+1) 



K = e ikn * I G(x;k)e lknx dx, (B18) 

-d- 



L = e _ ikn a i G{x;k)e- iknx dx, (B19) 



M = e -,fcn 3 / G{x]k)e lknx dx, (B20) 



tf = e ikn$ J 2 G(x;k)e- lknx dx, (B21) 



w 2n(l-r 2 e 42 * :a )' (n + 1) (1 - r 2 e * 2fe " d ) v ' 

The Langevin contribution is evaluated in the five regions. Since K, L, M and N depend linearly on the Fourier 
transform of the stochastic force operator, it should be noted that the coefficients also depend on the same way. 
In fact, since the stochastic force operator depend linearly on the bath's annihilation and creation operators, the 
Langevin contribution depend in that way too. 
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